The effect of long range interactions on the stability of classical and quantum solids 
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We generalise the celebrated Peierls' argument to study the stability of a long-range interacting 
classical solid. Long-range interaction implies that all the atomic oscillators are coupled to each 
other via a harmonic potential, though the coupling strength decays as a power-law l/x a , where x 
is the distance between the oscillators. We show that for the range parameter a < 2, the long-range 
interaction dominates and the one-dimensional system retains a crystalline order even at a finite 
temperature whereas for a > 2, the long-range crystalline order vanishes even at an infinitesimally 
^vq , small temperature. We also study the effect of quantum fluctuations on the melting behaviour of a 

one-dimensional solid at T — 0, extending Peierls' arguments to the case of quantum oscillators. 
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I. INTRODUCTION 

The role of the range of interactions in phase transitions and critical phenomena has been investigated for many 
years. For example, there are rigorous theorems which rule out the possibility of long-range ordered phases of a 
system with only short-range interactions [H . It is also well known that mean-field theory becomes exact in the limit 
of infinite range interactions (2j. Power-law interactions, which fall as l/x a with the increasing separation x between 
the interacting elements, can smoothly interpolate between these two extremes [1, H|. The implications of power-law 
interactions have been explored for several physical systems; these include, for example, classical Q and quantum 
spin glasses @, ionic systems jj], complex networks ifa], etc. Very recently, non-equilibrium phase transitions with 
long-range interactions have also received attention [9(. 

Power-law interactions can lead to non-trivial and counterintuitive results. For example, for the d-dimensional Ising 
[ model with ferromagnetic power-law interactions, the critical exponents depend on the range parameter a and the 
conventional short-range critical exponents are recovered for relatively large a (typically for a — > 2) 0, Hcf . The 
ferromagnetic to paramagnetic transition in a <i-dimensional Ising model with inverse-square interaction 1/x 2 is a 
O . special type of Kosterlitz-Thouless transition [ll[ driven by topological defects [13] . 

For many systems, heuristic reasoning, which are essentially extensions of Peierls' original arguments [l3(, have 
, been remarkably successful in predicting the (im-) possibility of long-range order. Some of these results were later 
^ ' established using rigorous analytical and numerical methods. Examples include predictions of the lower-critical 
ON ', dimension of the Ising model and the impossibility of the existence of a two-dimensional crystaline order at a finite 
temperature. 

In this work, employing Peierls'-type arguments, we examine the stability of long-range crystalline order in a d- 
. . dimensional system with power-law interactions. In other words, we consider a crystal at a finite temperature where 
the atomic oscillators interact via a harmonic potential with a strength that decays algebraically, i.e., as l/x a , with x 
being the seperation between the oscillators. It is to be noted that, the power-law form of interaction does not have 
an inherent length scale associated with it. This is in sharp contrast to the case where the strength of the interaction 
between two oscillators separated by a distance x decays exponentially as e~ x l L °. Clearly, in the latter case, L sets 
in a length scale over which the interactions are effective and therefore the melting behaviour turns out to be identical 
to oscillators with only nearest neighbour interactions. 

The Peierls' argument generalised to the power-law interactions suggests that the crystalline order persists for a one- 
dimensional solid even at a finite temperature, provided the range parameter a < 2 whereas for a > 2 the crystalline 
order is destroyed by infinitesimal thermal fluctuations. The case a = 2 corresponds to a marginal situation. We also 
comment on the stability of the crystalline structure at zero-temperature when the oscillators are quantum mechanical 
so that quantum fluctuations tend to destroy the crystalline order. 
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FIG. 1: (a) A One dimensional crystal and a (b) Two dimensional crystal. 
II. LONG RANGE INTERACTIONS 



Let us consider a one dimensional crystalline arrangement of harmonic oscillators with long-range interactions in 
the absence of any thermal or quantum fluctuations. The potential energy for this system(FigfT]) in the harmonic 
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approximation can be expressed as [141 ] , 



U = lA^^ina) q([n + m]a)Y (1) 



n m 



where K m is the effective spring constant, a is the interatomic seperation and n is the index representing the discrete 
sites on the crystai. The equation of motion for the atom at the position l na' is given by 

Mq\na) — K m [q([n + m]a) + q([n — m]a) — 2q{na)], (2) 

m>0 

where M is the mass of an oscillator. We assume the usual plane wave solutions of the type q(na) oc e l ( kna ~ wt ) ; to 
obtain the dispersion relation of the many-oscillator system given as 



Y,m>o K ^sin 2 (mka/2) 



^ 2 V M ® 

In the long wavelength limit, i.e. k—>0, which is relevant in determining the phase transition behaviour, Eqn.((3|) 
reduces to, 



E m >o^m 2 
M 



k\, (4) 



provided X)m>o m2 -^m converges. Now in the special case of nearest neighbour short range interactions, we recover 
the standard phononic 15] dispersion relation, to — ck, where c is the speed of sound. In the case of long range 

interactions, on the other hand, the interaction strength decays with distance as K m = — — (a >0), then holds 
for a > 4 since the sum in Eqn.([4]) converges. For 1 < a < 3, we get (see Appendix I) 

uj^ a -^l\ (5) 

It can be shown that the above dispersion relation holds for higher dimensional cases. We shall restrict our discussion 
to the range 1 < a < 3 to estimate the marginal dimension or the range of the stability of a crystalline structure. 

We shall now examine the effect of thermal fluctuations on the stability of the above solid using the formalism 
developed by Peierls' [l3j]. In the real space, the displacement of the atom at the position l na' can be written as a 
superposition of the different Fourier modes as, 

q(na) = J^e*™ (6) 

< k < 2ir/a, where the time dependence has been ignored. Making use of the classical equipartition theorem, we 
get 



ksT 

(kfcl 2 ) = (<?fc q-k) = -^—2, (7) 



The real space displacement, q(na) is immediately found to be, 



k R T r^kna 



(q 2 (na)) = J (q k q^ k )e lkna d d k = ^_ J ^~d% (8) 

where d denotes the dimensionality of space. Then (J 2 ) = (q 2 (na)) gives a measure of the fluctuation in the position 
of the atom at 'no! . To probe the marginality in the melting behaviour, we shall focus on the fluctuations in the 
displacement of the atom at the origin (na=0). 

Let us first enlist the form of (S 2 ) for solids with short range interactions in different spatial dimensions 

^lix/adk 



'JT f ford = 1 ' 
^/M ford = 2; (g) 

c-2-k I a k^dk 
J io -T2- ford = 3, 
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From the above Eqn.©, we see that (<5 2 ) diverges for a one and two dimensional solid in the limit k — > 0. Moreover, the 
2-d solid shows a log divergence, which means that the divergence is extremely slow. This log-divergence immediately 
establishes the spatial dimension d — 2 as the marginal dimension. For d = 3, on the other hand, (<5 2 ) is finite in 
the limit k — > 0, pointing to the existence of the crystalline order even at a finite temperature in a three-dimensional 
solid. 

Let us now extend the above analysis to solids with long range interactions. The form of the fluctation (S 2 ) is 
enlisted below: 



2-K/a 



(s 2 ) 



4 

J 



(■271 

J 



J J 



1< a < 3 



dk 


1 


k*- 1 


J, a -2 


dk 


1 


hot— 2 


k a ~ 3 


dk 


1 


k a ~ 3 


k a-4 



for d = 1, 



for d = 2, 



for d = 3, 



(10) 



The case with d = 1 turns out to be the most interesting as we observe that for all values a < 2, the crystal remains 
stable even at a finite temperature. At a = 2, we observe the first signature of a divergence. Therefore, we conclude 
that a — 2 happens to be the marginal case, with fluctuations diverging logarithmically with the system size. This 
settles a = 2 as the marginal case for a one-dimensional solid. If the long range interaction decays faster than 1/x 2 , 
the interaction is feeble and the solid fails to retain the crystal structure even at an infinitesimally small temperature. 



III. T=0 : EFFECTS OF QUANTUM FLUCTUATIONS 



In Eqii.©, if we put temperature T = 0, the average fluctuation (<5 2 ) vanishes. We now have to take care of the 
quantum fluctuations if the harmonic oscillators are quantum mechanical. In that case, we should now set the average 
energy of a mode to be the zero-point energy hwk of the corresponding oscillator. We obtain 



{ik q-k) 



Mw l 



Mwk ' 
(q k q- k )e lkna d d k 



M 



Akna 



-d d k. 



(11) 



(12) 



The above relations lead to a non-trivial modification in the marginal dimension of a solid with quantum mechanical 
oscillators, as illustrated below. In a spirit similar to that of the previous section, let us now examine the (S 2 ) for a 
quantum solid with short range interactions, 

r2n/adk 

J Jo Y fiwd =1, 



(s 2 



r2iv/a kdk 

Jj ° ~k~ 



for d 



(13) 



r 2tt /a k^dk 
"Jo -1- ioid 



Interestingly, under the effect of quantum fluctuations, both two and three dimensional crystals sustain the long 
range crystalline order, whereas the one dimensional solid shows a logarithmic divergence in the limit k — > 0. Hence 
in the case of quantum oscillators at T = 0, d = 1 is the marginal dimension whereas in higher dimensions a perfect 
crystalline order exists at absolute zero. The above finding when contrasted with the melting behaviour of a classical 
solid with short-range interacting oscillators, an interesting d quan t U m — *• (d + l) c iassicai correspondence, with respect 
to the marginal dimensionality arises. 
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TABLE I: Stability of Solids 





Short ran 


ge interactions 


Long Range Interactions 


Spatial 


Classical 


Quantum 


Classical 


Quantum 


dimension 


(T/0) 


(T=0) 


(T¥0) 


(T=0) 


d=l 


Unstable 


Unstable 


Stable(Ka<2) 


Stable(l<«<3) 






(In) 


Unstable(a>2) 


Unstable(ln) 










for (a>3) 


d=2 


Unstable 


Stable 


Stable(Ka<3) 


Stable 




(In) 




Unstable(ln) 










for q>3 




d=3 


Stable 


stable 


Stable 


Stable 



Let us now assume that the quantum oscillators are also long-range interacting as the classical oscillators of the 
previous section. A similar line of arguments leads to 



' r2ir/a dk 

~ J J 



1 



(■271 

J J 



fc(a-l)/2 


k( a 


-3)/2 


dk 




1 


k (a-3)/2 r 


k( a 


-5)/2 


dk 




1 


fc(a-5)/2 




-7)/2 



for d = 1, 



for d = 2, 



for d = 3, 



(14) 



1 < a < 3 



It is obvious from the above expressions that the question of any instability does not arise even in d = 1 if a < 3. 
Therefore a one dimensional long-range quantum solid retains the crystalline order for a finite strength of quantum 
fluctuations, if the decay of the long-range interactions with the distance between the oscillators is sufficently slow. 
When a = 3, 




which corresponds to the first sign of a divergence as a — > 3, and hence the situation is inferred to be marginal. The 
strength of the divergence keeps increasing till a = 4, when the (5 2 )~lnfc. 

We have summarised the different stability criteria obtained till now in Table I. 



IV. CONCLUSION 



In this work, we have analyzed the stability of classical and quantum solids using Peierls' argument, with a special 
emphasis on the one dimensional situation. We just note in the passing that mercury chain salts, in which mercury 
is intercalated into linear chains in AsF§ to give a compound Hg3_sAsFQ (S ^0.18 at 300K), is a system that closely 
resembles a one dimensional solid [l6| . 

In the presence of short range interactions, we never expect the 1-d solid to exist at any finite temperature (Eqn[5]). 
An infinitesimal thermal disturbance is sufficient to destroy long range order. The existence of long range interactions 
among the atomic oscillators, where the interaction strength falls off with distance as l/x a with a < 2, has a non- 
trivial effect on the stability of the one-dimensional crystal. Using a variant of Peierls' argument, we have shown 
that even a one-dimensional solid can sustain the crystalline order at a finite temperature if the range-parameter 
a is sufficiently small (a < 2). For a — 2, (S 2 ) shows a logarithmic divergence and hence this is the marginal 
case of stability. This is similar to what is observed also in the case of a ferromagnetic power-law interacting Ising 
model [l2l].That the inverse-square interaction is the marginal dimension even for melting of a long-range solid is an 
interesting observation. It should be mentioned that the situation a = 2 turns out to be a marginal situation in 
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various other long range interacting systems, for example in fracture models with variable range interactions |17l ] or 
in networks with long range links [181 ]. 

We have also explored the marginal dimension in the case where the oscillators are quantum mechanical. Using a 
generalised version of Peierls' argument extended to the quantum case, which to our knowledge is new, we observe 
an interesting d quan tum — * (d + I) classical correspondence (Table [XJ) . A two-dimensional short-range classical solid is 
marginally stable against thermal fluctuations whereas a one-dimensional short range quantum solid happens to be 
marginally stable against quantum fluctuations at T — 0. This reminds us of the well known correspondence between 
the critical behaviour of a d-dimensional quantum Ising model and [d + l)-dimensional classical Ising model [l9j . 
However, at an infinitesimally small temperature a one-dimensional quantum solid can not retain a long-range order, 
which means that the quantum effects are irrelevant at a finite temperature as the zero-point energy is vanishingly 
small compared to the thermal energy ksT . A unique correspondence emerges out of our study. We have shown 
that for a long-range one-dimensional classical solid, a — 2 is the marginal range of interaction. In the corresponding 
quantum case, a = 3 turns out to be the marginal case. There seems to exist a quantum-classical correspondence in 
the marginal range ofinteraction, just as in the case of dimension. 

One may also wonder about the existence of a tunable parameter in a quantum solid which may be varied to change 
the strength of the quantum fluctuations. Eqn.([7]), suggests that the mass of the atomic oscillator is an appropriate 
candidate for the above as the more massive the atom, the less are the fluctuations or uncertainty in its position. 

We should conclude with the comment that in deriving the marginal dimensions or ranges in above mentioned 
solids, we have used the Peierls' argument and its generalisation to the quantum case. Although, Peierls argument is 
fairly successful in predicting the marginal dimension or range in several other situations, it fails to throw any light 
on the transition temperature (for dimensions above the marginal dimension), nature of transitions or role of defects 
in the transition. However as we have shown, it predicts interesting results even in the case of melting of classical and 
quantum solids. 
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VI. APPENDIX I 

When the sum in Eqn.JU) is divergent, we recast Eqn.((3]) in the form of an integral as, 



2 l f™K m sin 2 {l3m)dm ^ 



where (3=ka/2. Now the integral is evaluated as follows(l < a < 3), 



'sin 2 ((3m)dm 1 [°° sin 2 (y)dy 

ZZ7, =P / 77, U 8 ) 



Since the second integral is just a constant, we obtain Eqn. ([!)]). For the special case a=3, the integrand in Eqn.JT 
diverges at m = 0. In this case, the dipersion relation turns out to be [14]], UJ= k\J\ In k\. 
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